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We present preliminary results for light, strange and charmed pseudoscalar meson physics from 
simulations using four flavors of dynamical quarks with the highly improved staggered quark 
(HISQ) action. These simulations include lattice spacings ranging from 0.15 to 0.06 fm, and sea- 
quark masses both above and at their physical value. The major results are charm meson decay 
constants /o, fj)^ and fojfo and ratios of quark masses. This talk will focus on our procedures 
for finding the decay constants on each ensemble, the continuum extrapolation, and estimates of 
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1. Introduction 

Lattice calculations of pseudoscalar meson decay constants, when combined with measured 
nonleptonic decay rates, can be used to extract CKM matrix elements and test the Standard Model. 
The lattice calculation includes a determination of the quark masses — fundamental parameters of 
the Standard Model. Here we present preliminary results for the charm meson decay constants /d 
and fo,, and the ratios of quark masses rric/ms and niu/nid. 

These calculations use lattices generated with the one-loop and tadpole improved Symanzik 
gauge action [1] and the HISQ fermion action [2]. We include four flavors of dynamical sea quarks: 
degenerate up and down, strange, and chaim. Iterated smearing in the HISQ action reduces taste 
violations by a factor of three relative to the asqtad action. An improved chaim-quark dispersion 
relation allows us to treat it the same way as the lighter quarks. Our lattice ensembles include 
ensembles with approximately physical light-quark masses, with volumes as large as 5.6^ fm^. 
For more details, see Refs [3, 4, 5]. The parameters of the ensembles used in this calculation are 
tabulated in J. Kim's talk [6]. In the first stage of analysis we fit pseudoscalar meson two-point 
con^elators to extract masses and amplitudes from random wall operators for a set of valence-quark 
masses. (See Ref. [6].) In the second stage of the analysis, we interpolate or extrapolate in valence- 
quark masses to find the tuned valence masses on each ensemble, and the decay constants evaluated 
at these valence masses, but at the sea-quark masses and lattice spacing of the individual ensemble. 
Then we combine the results of the different ensembles to make a continuum extrapolation and 
corrections for mistuned sea-quark masses. 

We use the pion decay constant fji to set the scale. Since this is determined from the same set 
of correlators as the charm meson decay constants, the analysis is self-contained in the sense that 
we do not need an intermediate quantity such as or r\ found in a separate analysis. This reduces 
the error from determining the lattice spacing, and simplifies the error analysis. Ensembles with 
physical light-quark masses dominate the analysis and allow a simple interpolation to adjust for 
light-quark mass corrections. However, it is likely that use of staggered chiral perturbation theory 
will lead to better controlled fits, and J. Komijani's talk describes progress in this direction [7]. 

2. Decay constants on each ensemble 

In the first stage of the analysis [6] we determine the masses and amplitudes for the two-point 
pseudoscalar correlators for a set of valence quark masses on each ensemble. These masses include 
0.9 and 1.0 times the sea charm quark mass, 0.8 and 1.0 times the sea strange-quark mass, and a 
range of lighter masses. Since the sea-quark masses were estimated before the production runs 
began, they are inevitably slightly mistuned. In the second stage of the analysis we determine 
tuned quark masses for each ensemble, and evaluate the decay constants at these valence-quark 
masses. Our primary analysis uses fji to fix the lattice spacing, and proceeds as follows. To 
estimate statistical errors, this whole procedure is done inside a jackknife resampling. 

1. Interpolate or extrapolate in the light valence-quai^k mass m/ to the point where M^/ fji has 
its physical' value. This fixes ami and the lattice spacing a. 

2. Interpolate or extrapolate in valence strange-quark mass arris to where 2M|- — has its 
physical value. This fixes rris. 

'adjusted for E&M and finite size — see later 
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Figure 1: fo and fo^ on the different ensembles. Fits that are quadratic and linear in a^, or use only physical 
quark mass ensembles are superimposed on the data. More details are given in the text. 

3. Interpolate or extrapolate to the charm valence mass where Md, is correct. This fixes ruc. 

4. Find rtid — niu from E&M adjusted —K^ mass difference [8] 

5. Find, by interpolation or extrapolation, /k at the adjusted up-quark mass. 

6. Find /o and (a check) at the adjusted down and charm masses. 

7. Find /o^ at the adjusted strange and charm masses. 

8. Find M^^. (check) at the adjusted charm-quark mass. 

We have also used a tuning procedure where we replaced fj^ and Mj^ by the decay constant, 
fp4s, and mass, Mp4s, of a fictitious meson with degenerate valence quark masses at 0.4 times 
the strange-quark mass, using the values for fp4s and Mp4s in MeV determined from three flavor 
asqtad data. Finally, we varied the procedure by first fixing the lattice spacing from the static-quark 
potential, and then tuning the valence quark masses as described above. 

In Eq. 2. 1 we list results from the above procedure for the most influential ensemble (a 0.09 
fm with physical masses). These results are adjusted for finite size and and some electromagnetic 
effects, but are not adjusted for mistuned sea-quark masses or extrapolated to the continuum limit. 
The errors here are statistical only but, as discussed above, errors from scale setting and valence- 
quark mass tuning are included in the statistical error here. The D+, D° and Tj^ masses are not used 
in the tuning, so they can be compared to their experimental values, shown in parentheses. 

a = 0.08792(10) fm am, = 0.001333(5) am, = 0.03648(11) am, = 0.4323(7) 
m„/mrf = 0.480(6) =27.36(3) mc/m, = 11.851(17) 

fK= 155.12(22) MeV 

Mdo = 1867.8(1.5) MeV (cf 1864.8) M^^ = 1870.2(1.2) MeV (cf 1869.6) 

M^^ = 2980.2(3) MeV (cf 2980.3(1.2)) 

/d = 208.75(1.03) MeV /d, = 246.60(20) MeV /^y/^ = 1.181(6) i) 

3. Continuum extrapolation 

The results for fjj and fo^ on each ensemble are shown in Fig. 1 . In this figure the physical 
quark mass ensembles, labelled as "m/ = 0.04ms", have the smallest statistical errors. This is 
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Figure 2: The ratios foj/o (left panel) and mc/m^ (right panel) on the diftferent ensembles. The superim- 
posed fit functions are the same forms as in Fig. 1. 

due to their larger physical volumes, where the random wall and Coulomb wall sources, together 
with the average over spatial position of the sink operator, result in better statistics. EiTors ai^e 
large on some of the m/ = 0.2 wis ensembles because the lightest valence-quark mass used on these 
ensembles was 0.1 m^., so a significant extrapolation in valence light-quark mass was necessary. 
When these points are fit to functions of lattice spacing and sea-quark mass, the physical quark 
mass points dominate the fit. Indeed, we could get reasonable results by simply using the physical 
quark mass ensembles. In practice, we use the physical quark mass and the mj =0.1 ensembles, 
where the mi = 0.1 nig ensembles allow correction for mistuning of the light sea-quark masses in 
the physical quark mass ensembles and, because of the m/ = 0.1 ensemble at lattice spacing 
a « 0.06 fm, help to determine the dependence on the lattice spacing. 

Figure 1 shows three different fits to this data. The magenta line uses all four of the lattice 
spacings, and is quadratic in and linear in Ught sea-quark mass. The straight cyan line in Fig. 1 
omits the a = 0.15 fm points, and is linear in both and ni/ sea- The green line is a quadratic fit 
using only the physical quark mass points, with small adjustments to correct for sea-quark mass 
mistuning. The symbols at a = are the continuum limits of these fits, showing the statistical 
error. The size of these symbols is proportional to the p-value of the fit, with the symbol size in 
the graph legend corresponding to 50%. Since we expect both order and corrections to the 
data, and since some other accurately determined quantities such as the mass of the rjc clearly need 
terms to fit the results, we take the quadratic fit results as our central value, with the larger of 
the difference between the quadratic and linear- fits or the difference between the quadratic and 
physical-quai^k-mass-only fits as an estimate of the systematic eiTor coming from our choice of fit 
form. 

Figure 2 shows, with the same notation, the ratios foJ/o and nic/ms on each ensemble, to- 
gether with fits to the same functional forms and ranges of data as in Fig. 1 . 



4. Systematic errors 

Statistical errors from the jackknife analysis incorporate errors in lattice spacing determina- 
tion and valence-quark mass tuning, since this tuning is redone for each jackknife sample. The 
following additional systematic errors are included in our eiTor budget. 
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Figure 3: Comparison of lattice results for fjj and fo^ (left panel) and for fj)^ / fo (right panel). Results 
are from Ref. [10] and this work. Diamonds represent results with 2, octagons with 2+1, and squares with 
2+1+1 dynamical flavors. For this work the red error bars (lower) are statistical only, while the blue error 
bars (higher) include systematic error estimates. 

Effects of excited states in the two-point coirelators were estimated by varying the time range 
in the fits and the priors for excited state masses with reasonable ranges [6]. 

Effects of finite spatial size were tested by running otherwise identical ensembles with three 
different spatial sizes. This was done at a 0. 12 fm with nij/ms = 1/10. The results, summarized 
in Fig. 4, show that the significant effects are in the light pseudoscalar sector, on Mj^, and fx, and 
that chiral perturbation theory is a reasonable guide. From these results, we estimated the effects 
on Mji, fji and at the physical light-quark mass in a 5.5 fm box. (Since our fits are evaluated at 
this mass, it doesn't matter that the ensembles with other sea-quark masses would have different 
finite-size effects.) Then, in the analysis procedure described in section 2, we used the values of 
Mji, fji and /k adjusted to this size box. (That is, we compute fx/fniS-S fm) etc.) Afterwards, we 
use these same factors to coiTcct the computed /k back to its infinite volume value. We use the 
values with the finite-size adjustments as our central values, and take one half of the shifts when 
these adjustments are omitted as a remaining systematic error. 

The effects of electromagnetic interactions and isospin breaking (m„ / m^/) are considered to- 
gether. We tune using the 71° mass, where electromagnetic and isospin breaking effects are expected 
to be small. Then we use a separate calculation of electromagnetic effects on the kaon mass [8], 
whose main result is parameterized in terms of a violation of Dashen's theorem 
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Here M/f + adj is the mass adjusted to remove the effects of electromagnetism, and the result of 
Ref. [8] is Aem = 0.65(26). Then, in tuning the strange-quark mass from 2M|- — M^, we use the 
average squared kaon mass M\ = {m^+ +M^o^ /2. We also use this result combined with our 
meson mass measurements to determine the up-down quark mass difference. 
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Figure 4: Spatial size effects on M^i, Mk, fn and /k (left side), and on Mo, Mo,, fo and fo, (right side). 
Solid lines are the one loop chiral perturbation theory forms. To show the magnitude of the effects, green 
error bars show an arbitrary value ±1 MeV, and magenta error bars ±1%. 



Source 


fo 


fo, 


foj fo 




rriu/md 


Statistics 


3.0 MeV 


0.5 MeV 


0.016 


0.041 


0.009 


Excited 


0.5 MeV 


0.5 MeV 


0.004 


0.005 


0.003 


Volume 


0.3 MeV 


0.2 MeV 


0.0004 


0.008 


0.0005 


E&M 


0.04 MeV 


0.12 MeV 


0.001 


0.001 


0.021 


Scale setting 


2.0 MeV 


1.3 MeV 


0.003 


0.02 


0.023 


fit form 


2.9 MeV 


3.3 MeV 


0.010 


0.09 


0.011 



Table 1: Error budget: The statistical error includes error from valence-quark mass tuning and from a 
continuum extrapolation with fixed fit form. Estimation of the systematic error is discussed in the text. 



where is obtained from the difference in between the two lightest valence-quark masses, 
with the strange mass set at the sea strange mass. We estimate the remaining systematic errors 
from electromagnetism by combining, in quadrature, the effects of changing Aem by one standard 
deviation with the effects of shifting the average kaon mass squared by 900 MeV^, which is the 
full result of a preliminary computation of electromagnetic effects on this quantity [8]. When 
computing fo, the Ught valence-quark mass is interpolated to the resulting down-quark mass, and 
similarly for other quantities involving valence up and down quarks. We expect that the effects of 
electromagnetic interactions and isospin violation in the sea quarks are small, and ignore them. 

Errors from the scale setting and quark mass tuning procedure are estimated from the change 
in results when ri or fp4s are used to set the scale instead of fji. 



5. Conclusions 

We find preliminary results fo = 209.2(3.0) (3.6) MeV, fo, = 246.4(0.5) (3.6) MeV, /d,//d = 
1.175(16)(11), m,./m, = 11.63(4)(9) and m„/md = 0.505(9)(33). Our results are consistent with 
those of other collaborations [10], and the errors are comparable to the most precise published 
calculation by HPQCD. In the future we will further reduce the eiTor by using staggered chiral 
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perturbation theory, which will allow us to make better use of the correlators with unphysical 
valence-quark masses [7] , and by the addition of an ensemble with physical sea-quark masses at a 
lattice spacing of 0.06 fm. 
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